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1.0  INTRODUCTION 

Current  techniques  for  the  prediction  of  propulsion  system  drag  for  aircraft  rely  heavily 
on  wind  tunnel  tests.  Yet,  these  tests  do  not  incorporate  complete  simulation  of  the  jet 
exhaust.  Preliminary  configuration  tests  in  many  cases  use  solid  plume  simulators  (Ref.  1). 
Validation  of  the  final  design  is  accomplished  in  the  wind  tunnel  with  high-pressure  air  to 
simulate  the  nozzle  exhaust  flow.  However,  viscous  interaction  effects  for  a  cold  jet  differ 
from  those  for  a  hot  jet  since  mixing  between  the  nozzle  exhaust  and  local  external  flow 
depends  largely  on  temperature,  species,  and  velocity  gradients.  Experimental  and  analytical 
investigations  at  the  Arnold  Engineering  Development  Center  (AEDC)  (Refs.  2  through  6) 
have  attempted  to  characterize  the  various  simulation  parameters  necessary  to  obtain  valid 
test  results,  but  it  is  still  not  possible  to  accurately  predict  the  afterbody  drag  on  nozzle- 
afterbody  models.  If  a  sufficiently  powerful  analytical  or  numerical  tool  were  available,  the 
proper  procedure  for  simulating  a  hot-jet  exhaust  flow  with  a  cold  fluid  might  be  developed. 
Creation  of  this  tool  was  the  objective  of  the  present  investigation. 

Current  methods  (Refs.  6  through  14,  for  example)  for  predicting  nozzle-afterbody  drag 
have  relied  on  patched  viscous-inviscid  solutions  wherein  the  complete  flow  field  is 
computed  in  an  iterative  fashion,  alternating  between  a  potential  flow  solution  and  a  viscous 
solution.  With  one  notable  exception  (Ref.  14),  these  methods  require  severe  under¬ 
relaxation  of  the  viscous  displacment  correction,  and  the  solutions  obtained  are 
unfortunately  highly  dependent  on  the  relaxation  technique.  In  spite  of  this  ditticulty,  some 
of  the  methods  are  in  good  agreement  with  experimental  data,  particularly  the  method  of 
Presz,  King,  and  Bureau  (Ref.  13). 

Another  approach  to  the  afterbody  problem  is  that  of  Holst  (Ref.  15)  who  developed  a 
numerical  solution  of  the  time-dependent,  compressible  Navier-Stokes  equations.  Holst 
demonstrated  excellent  agreement  with  surface  pressure  data  obtained  from  experiments 
with  circular-arc  afterbodies  and  solid  cylindrical  plume  simulators. 

The  present  study  is  also  based  on  numerical  solution  of  the  Navier-Stokes  equations. 
The  planar  two-dimensional  computer  program  developed  by  Deiwert  (Ref.  16)  was 
modified  to  enable  computation  of  axisymmetric  flow  over  nozzle  afterbodies  with  either 
real  or  simulated  plumes.  The  basic  numerical  algorithm  is  the  explicit  predictor-corrector 
scheme  of  MacCormack  (Ref.  17).  Turbulence  modeling  is  accomplished  using  the 
algebraic,  two-layer  eddy  viscosity  model  of  Baldwin  and  Lomax  (Ref.  18)  with  a  novel  local 
modification  as  a  function  of  vorticity.  The  results  presented  are  limited  to  comparisons  of 
computations  with  the  recent  detailed  flow-field  measurements  of  Benek  (Ref.  19)  and  with 
the  surface  data  of  Reubush  (Ref.  1). 
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2.0  NUMERICAL  METHOD 

2.1  AXISY  MMETRIC  COMPRESSIBLE  FLOW  EQUATIONS 

The  equations  describing  the  axisymmetric  flow  of  a  compressible,  viscous  fluid  without 
body  forces  can  be  written  in  weak  conservation  form  as  follows  (Refs.  20  and  21): 
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The  density,  g,  is  a  mean  value  and  the  velocities,  u  and  v,  are  mass-averaged  values.  The 
mass-averaged  specific  integral  energy,  e,  is  related  to  the  mean  total  energy  per  unit  volume, 
E,  by 


(2) 


and  the  perfect  gas  equation  of  state 


P  =  (y  —  ])  pe  (3) 

is  used  to  define  the  mean  pressure,  P.  The  turbulence  model  (see  Section  2.2)  expresses  the 
Reynolds  stress  tensor  in  terms  of  the  scalar  eddy  viscosity,  e,  which  is  additive  to  the 
molecular  viscosity,  g,  but  not  incorporated  in  the  second  coefficient  of  viscosity,  X,  which  is 
defined  as  X  =  -2/3  g.  The  Sutherland  viscosity  relation  is  used  and  the  ratio  of  specific 
heats,  y,  is  maintained  constant  at  1 .4.  Laminar  and  turbulent  Prandtl  numbers  are  lixed  at 
Pr  =  0.72  and  PrT  =  0.9,  respectively. 


2.2  TURBULENCE  MODEL 

The  two-layer  algebraic  turbulence  model  of  Baldwin  and  Lomax  (Ref.  18)  was  selected 
for  this  study  because  its  use  eliminates  the  necessity  of  locating  the  boundary-layer  edge. 
The  required  length  scale  is  determined  from  the  distribution  of  vorticity  in  the  vicinity  of 
the  body  surface.  The  eddy  viscosity,  e,  is  defined  as  a  function  of  the  normal  distance  from 
the  wall,  y,  in  the  inner  region  as 


and  in  the  outer  region  as 


fout  =  kccpp  fwake  fkleb  ^ 

The  inner  form  is  used  out  to  the  point  where  ein  =  eom  and  from  that  point  e 
auxiliary  functions  appearing  in  Eqs.  (4)  and  (5)  are  defined  as  follows: 
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The  constants  appearing  in  this  formulation  are  the  following  values  assigned  by  Baldwin 
and  Lomax: 


A+  =  26 
CCP  =  1.6 

CKLEB  =  03 
CWK  =  0.25 

k  =  0.4 
K  =  0.0168 


Holst  (Ref.  15),  Shang  and  Hankey  (Ref.  22),  and  others  have  shown  that  turbulence 

^  - -  — - - - ~ 

models  derived  from  turbulent  boundary-layer  data  at  equilibrium  conditions  do  not  yield 

good  results  when  applied  to  separated  flows.  Their  approach  to  introducing 

nonequilibrium  effects  wras  to  arbitrarily  _select  a  streamwise  location  as  a  reference 
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condition  and  exponentially  relax  the  calculated  eddy  viscosity  through  the  separated  region. 
This  approach  was  included  in  the  present  analysis  but  yielded  disappointing  results. 

One  characteristic  of  boundary  layers  in  equilibrium  (observed  in  the  present  study)  is 
that  the  streamwise  gradient  of  vorticity  is  small.  It  follows,  therefore,  that  the  degree  of 
nonequilibrium  effects  can  be  quantified  by  examining  the  vorticity  gradient.  This  concept  is 
used  to  locally  adjust  the  Baldwin-Lomax  eddy  viscosity  in  the  form 

h  AX  dl(t)l  l 

'»  ■  '  v  ~I"T  d,  )  (13) 

where  5  is  an  arbitrary  constant  of  order  unity.  As  a  matter  of  choice,  the  maximum  change 
in  e  is  limited  to  ±  100  percent.  This  formulation,  termed  “vorticity  gradient  relaxation,”  is 
new  and  has  not  been  fully  explored. 

2.3  SOLUTION  ALGORITHM 

For  most  practical  situations  the  computational  mesh  is  nonorthogonal  and  it  is 
convenient  to  recast  Eq.  (1)  into  an  integral  form  (Ref.  23)  valid  for  each  volume  element  as 

-2-  f  rUdVol+  / rH-i^s=  f  A  dVol 

<?‘  vTl  5  Vol  (14) 

where  H  is  the  vector  sum  of  F  and  G  and  n  is  the  outward  unit  vector  normal  to  the  surface 
of  the  computational  cell.  The  area  vectors  nds  are  represented  as  S  as  shown  in  the  sketch  of 
an  arbitrary  quadrilateral  volume  element. 
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Early  versions  of  the  present  formulation  used  the  “splitting”  technique  of 
MacCormack  and  Paullay  (Ref.  23)  which  worked  quite  well  provided  the  mesh  was  aligned 
with  the  radial  coordinate.  Special  applications,  particularly  inviscid  solutions  with  irregular 
geometries,  indicated  that  an  unsplit  solution  procedure  would  be  more  stable  and  robust. 
The  splitting  technique  was  discarded,  in  spite  of  the  time-saving  attribute,  because  of 
pronounced  oscillations  in  the  dependent  variables  between  predictor  and  corrector  steps. 


The  MacCormack  algorithm  (Ref.  17)  as  applied  to  Eq.  (14)  may  be  written,  with 
superscripts  p  and  c  representing  predictor  and  corrector,  as  follows: 
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(15) 


The  forward  and  backward  permutations  of  the  indices  i.j  are  varied  cyclically  each  time 
step  (total  of  four  steps)  to  balance  the  truncation  error  at  each  cell.  Thus,  the  dependent 
variables  U  are  explicitly  advanced  in  time  throughout  the  mesh  to  Up,  thence  corrected  to 
UC 


The  centroid  radius,  rB,  is  currently  the  average  of  the  radii  at  the  cell  corners.  An  earlier 
version  of  the  program  utilized  the  correct  radial  weighting  to  define  rB,  but  yielded  no 
discernible  difference  relative  to  averaging. 

The  maximum  allowable  time  step.  At,  is  given  by  the  Courant-Friedrich-Lewy  (CFL) 
stability  criteria.  Ref.  21. 

2.4  FINE-MESH  CONSIDERATIONS 

The  CFL  stability  criteria  require  prohibitive  computation  time  in  regions  of  closely 
spaced  mesh  points  which  are  necessary  to  resolve  the  viscous  region  adjacent  to  a  surface. 
Again,  recourse  was  made  to  MacCormack’s  initiative  (Ref.  24)  and  Eq.  (1)  recast  into  time- 
split  hyperbolic  (inviscid)  and  parabolic  parts.  The  hyperbolic  equations  are  solved  using  an 
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explicit  numerical  method  based  on  characteristic  theory  and  the  parabolic  equations  solved 
implicitly  using  a  simple  tridiagonal  algorithm.  The  fine  mesh  is  constrained  to  be  stationary 
and  parallel  to  the  surface  and  Deiwert’s  (Ref.  25)  local  coordinate  rotation  incorporated  to 
improve  the  solution  accuracy. 


Equation  (1)  is  locally  transformed  from  (x,r)  to  the  surface-oriented  (£,y)  coordinates  to 


yield 
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at 
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Equation  (16)  is  then  split  with  the  ^-derivative  terms  handled  as  with  Eq.  (15)  but  the 
y-derivative  terms  further  split  into  viscous  and  inviscid  terms.  The  nomenclature  associated 
with  the  coordinate  rotation  is  defined  in  the  following  sketch 


f  -  £2 


and  tangential  and  normal  velocity  components  defined  as 
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The  hyperbolic  inviscid  equations  then  may  be  written  as 
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which  may  be  reduced  to  the  form 
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where  the  continuity  and  energy  equations  are  redundant.  For  small  v ' ,  the  last  expresssion 
of  Eq.  (19)  may  be  temporarily  neglected  to  obtain  characteristic  relations 

J  X  =  v'  ±  c 
dt 
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(20) 


where  the  speed  of  sound  is  given  by 


(21) 


The  normal  velocity,  v',  and  pressure,  P,  are  provisionally  advanced  in  time  using  the 
characteristic  relations  (evaluating  the  term  c2v/r  at  the  previous  time  value),  the  tangential 
velocity  is  updated  via  the  last  expression  of  Eq.  (19),  and  the  density  is  updated  using  dp  = 
dP/c2.  These  estimated  variables  are  then  used  to  calculate  the  inviscid  flux  vectors  Hj  •  S| 
and  H3  •  S3  of  Eq.  (15)  which  are  held  fixed  during  both  the  predictor  and  corrector  steps. 
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The  parabolic  viscous  equations  may  be  written  as 


where  the  continuity  equation  is  represented  by  stationary  density  (and  a*  - 
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dul  =  l  dxT\  s*°6 
(?t  rpJS  dy  rp 

jv_\  t  jrT2  , 
d<  rpJS  dy  rp 


where 


S(p  +  f)  du '  S(p  -f  f)  \  du  dr_  dv_ 


dy 


J  )  dy  dy  d£ 


,  9  /c'du  S'  dv  \  (c'  dj_  _  t -  daN  ( 

+  2^ar  *%){*&  "3f  dy  /  j 

Ti  A  -W  A' s  dv'-  .AS  (dr  du  dx  dv  \ 

T2  =  ^+2(p  +  ,)jT  —  +T^aT  H  dy 

2(a  +  t)S  A  ,  $v  c-d*\(<Z'd\i  .  S'dy\ 
- J—  (^dT  “  r  57;  ft  +  Sr  ) 


+  + 


(22) 

x  -  P,  etc.). 

(23) 
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These  equations  are  solved  numerically  (see  Ref.  24)  using  first-order  Laasonen  implicit 
differencing  and  tridiagonal  reduction.  The  energy  equation  is  replaced  by  two  kinetic 
energy  relations  and  an  internal-energy  equation  written  as 


Solution  of  Eq.  (24)  requires  only  one  additional  tridiagonal  reduction. 

The  solution  of  the  fine  mesh  is  advanced  in  time  using  basically  two  operators  in  the 
sequence  L!(At/2)L2(At)Li(At/2).  The  operator  L|  is  Eq.  (15)  without  the  viscous  terms  in 
the  “j”  direction  (unsplit)  and  the  operator  L2  is  the  implicit  split  viscous  solution,  Eqs.  (23) 
and  (24). 

2.5  SMOOTHING 

Two  types  of  smoothing  are  incorporated  to  eliminate  “wiggles.”  The  basic 
MacCormack  algorithm  is  subject  to  nonlinear  stability  problems  and  his  expansion¬ 
averaging  technique  (Ref.  26)  is  used  to  maintain  stability.  The  second  type  of  smoothing 
used  was  also  suggested  by  MacCormack  and  consists  of  explicitly  including  some  of  the 
third-order  truncation  errors  in  evaluating  the  flux  vectors  H  •  S  in  Eq.  (15). 

The  source  of  the  third-order  smoothing  terms  may  be  seen  by  considering  the  following 
model  nonlinear  equation 


dL  *  *£ 

dt  dr 


0 


(25) 


with 


£  =  f(x,t) 

g  =  gW 
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and  applying  the  MacCormack  algorithm  in  the  form 


f?  -  f"+1  =  %<l 


i+fi  ~t  [(fs)i"(f  s)i-i] 


The  solution  is  considered  exact  at  time  n  and  Taylor  series  expansions  in  the  form 


(fg)i+1  =  <fs>i  +  Ax^f- 

©i_i  = 


,3  J  3 , 


21  dx2  3! 

<?x3 

Ax2  Ax3 

& 

CO 

w 

X 

fM 

dx3 

(7.  -  f. 

)2  d2S 

(?  .  -  f . 

)s  a!* 

+  tt 

dt2 

1  3t 

at3 

(26) 


(27) 


are  used  to  evaluate  the  error  encountered  at  time  n  +  1  by  comparison  with  the  Taylor  series 
expansion  of  the  exact  solution.  Defining 

*  -  f  ir +  B  (28) 


the  MacCormack  algorithm  yields 


At3  d_ 
4  <5x 


At2Ax 

4 


ArAx  d  ( thfe  die  \  _  AtAx  ~ 
2  5i/  6 


and  the  exact  solution  is 


f  (t  +  At) 


c/x 


All  A  L 

2  dx  x  dx ) 


(29) 
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V.. 

Comparison  of  these  expressions  shows  that  MacCormack’s  algorithm  on  a  uniform  mesh  is 
second-order  accurate  for  this  class  of  nonlinear  equations.  Indeed,  much  of  the  “flavor”  of 
the  third-order  terms  is  present  and,  with  Ax  =  4>At,  the  truncation  error  at  third  order  is 
given  approximately  by 

fP+ 1  _  f  (t  +  At)  =  —  ( ^  2  +  M  _  AC  ,2  r?afg 

1  6  V*/  6  ax2  6  a*3  (3i) 

Since  these  terms  are  small  it  is  not  necessary  to  code  them  exactly  and  note  simply  that  the 
truncation  errors  are  proportional  to  the  first,  second,  and  third  spatial  derivatives.  In 
practice,  the  CFL  stability  criterion  for  the  split  Euler  equations  define  the  relationship 
between  At  and  Ax  which  specifies  4>  -  u  +  c  or  |v|  +  c. 

2.6  COMPUTATIONAL  MESH  AND  BOUNDARY  CONDITIONS 

Representative  mesh  constructions  are  given  in  Fig.  1.  A  fine  exponentially  stretched 
mesh  is  constrained  to  be  parallel  to  the  body  surface  and  extends  into  the  wake/plume 
region,  if  present,  as  shown  in  Fig.  lb  with  the  first  node  point  within  the  viscous  sublayer. 
Coarser  mesh  point  spacing  is  specified  in  regions  of  reduced  viscous  effects  with  radial 
exponential  stretching  in  the  outer  flow  and  proportionate  radial  spacing  within  the  plume. 
The  mesh  includes  points  at  the  boundaries  for  specification  of  boundary  conditions. 

Zero-slip  boundary  conditions  at  viscous  surfaces  are  applied  using  the  reflection 
principle  on  velocities,  evaluation  of  the  normal  momentum  equation  for  pressure,  adiabatic 
wall  conditions  for  internal  energy,  and  the  equation  of  state  for  density.  Inflow  boundary 
conditions  are  fixed  in  time  at  the  initial  condition  as  is  the  upper  boundary  condition.  The 
outflow  or  downstream  boundary  condition  is  updated  assuming  zero  gradients  in  the 
streamvvise  direction  and  the  centerline,  if  present,  is  treated  with  symmetry. 

Initial  conditions  consist  of  uniform  flow  everyw  here  in  the  coarse  mesh  and  a  boundary- 
layer-type  profile  in  the  fine  mesh  obtained  from  Whitfield’s  velocity  profile  (Ref.  27)  with 
specification  of  the  desired  inflow  boundary-layer  displacement  and  momentum  thicknesses. 
For  the  real-plume  calculations,  the  interior  inflow  is  computed  from  isentropic  relations 
with  stagnation  pressure  and  temperature  being  smoothly  varied  to  the  desired  conditions 
over  the  first  100  time  steps. 
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3.0  RESULTS  AND  DISCUSSION 

3.1  FLOW  OVER  A  BOATTAIL  WITH  SIMULATED  PLUME 

Benek  (Ref.  19)  has  recently  completed  an  experimental  investigation  specifically 
designed  to  provide  information  for  the  validation  of  axisymmetric-flow  computer 
programs.  The  geometry  used  in  that  program  is  represented  in  Fig.  la  and  consists  of  a 
circular-arc  boattail  (f/D  =  0.8)  followed  by  a  contoured  solid  plume  simulator.  Detailed 
flow-field  measurements  using  iaser  velocimetry  were  made  including  results  at  locations 
corresponding  to  the  boundary  of  the  computational  domain. 

Computed  results  from  the  present  analysis  are  compared  with  Benek’s  data  in  Figs.  2 
and  3.  The  surface  pressure  comparison  (Fig.  2)  show's  good  agreement  except  in  the  vicinity 
of  the  cusp  where  separated  flow  is  evident.  (If  the  measurements  were  to  be  adjusted  from 
the  stated  =  0.64  to  an  apparent  free-stream  Mach  number  of  0.635,  then  the 
computations  and  data  would  be  in  excellent  agreement  upstream  of  the  boattail.) 
Application  of  vorticity-gradient  relaxation  (5  =  2)  to  the  Baldwin-Lomax  turbulence  model 
(5  =  0)  significantly  changes  the  computed  results.  The  effect  of  relaxation  is  vividly  shown 
in  the  velocity  profile  comparisons  given  in  Fig.  3.  Relaxation  improves  the  solution  at  some 
locations  (at  separation  and  reattachment  in  particular)  but  has  adverse  effects  at  other 
stations.  It  is  obvious  that  the  turbulence  model  could  be  tuned  to  yield  computed  results  in 
complete  agreement  with  this  experimental  data,  and  equally  obvious  that  the  turbulence 
model  is  the  weakest  link  in  the  computation  of  turbulent  separated  flow. 

3.2  EFFECT  OF  THIRD-ORDER  SMOOTHING 

The  flow  Held  about  the  boattail  is  presented  in  Fig.  4  in  the  form  of  isobars  to  illustrate 
the  smoothing  achieved  by  explicit  consideration  of  the  third-order  truncation  errors.  The 
“wiggles”  evident  in  Fig.  4a  are  suppressed  in  Fig.  4b  by  including  the  gradient  smoothing 
term,  dfg/dx,  of  Eq.  (31).  Although  not  shown,  the  gradient  term  is  particularly  effective  in 
smoothing  a  solution  containing  embedded  shock  wraves.  Further  addition  of  the  curvature 
smoothing  term,  d2fg/3x2,  suppresses  the  wiggles  over  the  cusp  region  as  shown  in  Fig.  4c. 
Since  this  solution  is  sufficiently  smooth,  the  third  derivative  term  was  not  included  in  the 
computer  program.  A  note  of  caution  is  in  order  —  this  type  of  smoothing  is  so  powerful 
that  wiggles  arising  from  coding  errors  can  also  be  smoothed. 
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3.3  FLOW  OVER  A  BOATTAIL  WITH  REAL  PLUME 

Reubush  (Ref.  1)  presents  surface  data  for  a  variety  of  boattails  with  and  without 
cylindrical  plume  simulators.  One  configuration  (No.  1)  was  selected  for  comparison  with 
the  present  program  and  the  results  are  given  in  Fig.  5.  The  pressure  distribution  over  the 
boattail  is  reasonably  well  predicted  but  it  is  not  known  to  what  degree  of  accuracy  the 
plume  is  calculated.  The  computed  velocity  vectors  at  the  nozzle  exit  are  given  in  Fig.  6  to 
illustrate  the  flow-field  complexity.  It  is  highly  desirable  that  some  agency  obtain  detailed 
experimental  measurements  about  this  type  of  nozzle  afterbody  with  a  real  plume  to  provide 
code  verification  data.  Modeling  of  the  mixing  and  entrainment  effects  of  the  plume  cannot 
be  accomplished  by  comparison  with  surface  pressure  data  alone. 

4.0  CONCLUDING  REMARKS 

A  generalized  axisymmetric  version  of  Dei wert’s  compressible  Navier-Stokes  computer 
program  has  been  developed.  The  immediate  objective  was  the  generation  of  a  tool  for  the 
prediction  of  jet  simulation  parameters  to  be  used  in  wind  tunnel  tests  of  nozzle-afterbody 
configurations.  Limited  comparisons  of  computed  results  with  experimental  data  indicate 
that  the  perfection  of  such  a  tool  awaits  solution  of  the  turbulence  modeling  problem. 
Nonetheless,  for  a  given  turbulence  model  it  is  expected  that  the  computer  program  can  be 
used  to  reasonably  well  predict  the  incremental  effects  on  nozzle-afterbody  drag  as 
influenced  by  internal  geometry,  nozzle  pressure  ratio,  and  exhaust  temperature. 
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I  I 


b.  Real  plume 

Figure  1.  Representative  computational  meshes. 


VERTICAL  DISTANCE  FROM  SURFACE  (r-R)/D 


•  DATA  OF  BENEK  (REF.  19) 

-  BALDWIN  -  LOMAX  (REF.  18)  TURBULENCE  MODEL 

— —  VORTICITY -GRADIENT  RELAXATION  (»-2) 


LONGITUDINAL  VELOCITY,  o/u. 


Figure  3.  Comparison  of  computed  and  measured  velocity  profiles 
about  an  afterbody  with  simulated  plume. 
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a.  Expansion-averaging  only 


b.  Gradient  smoothing 


c.  Curvature  smoothing 
Figure  4.  Effect  of  third-order  smoothing. 
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O  DATA  OF  REUBUSH  (RER  I) 

NPR  =  i  I)  COMPUTED  (8  =  2) 
- N  PR  =  3.6  J 


Figure  5.  Comparison  of  computed  and  measured  surface  pressure 
coefficients  on  an  afterbody  with  real  plume,  M,„  =  0.8. 
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NOMENCLATURE 

Source  term  in  Eq.  (1),  also  van  Driest’s  constant 
Speed  of  sound,  also  model  chord 
Turbulence  model  Baldwin-Lomax  constant  (see  Eq.  (5)) 
Turbulence  model  Klebanoff  constant  (see  Eq.  (9)) 
Pressure  coefficient,  (P-P0o)/q* 

Turbulence  model  wake  constant  (see  Eq.  (10)) 

Model  diameter 

Mean  total  energy  per  unit  volume 
Specific  internal  energy 
See  Eq.  (1) 

Turbulence  model  (see  Eq.  (9)) 

Turbulence  model  (see  Eq.  (1 1)) 

Turbulence  model  (see  Eq.  (10)) 

General  functional 
See  Eq.  (1) 

General  functional 
Jacobian 

Turbulence  model  constant 
Turbulence  model  constant 

Turbulence  model  mixing  length,  also  model  length 

Mach  number 

Unit  normal  vector 

Static  pressure 

Prandtl  number,  laminar 
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PrT  Prandtl  number,  turbulent 

Q  Heat 

q  Dynamic  pressure 

R  Model  radius 

Re  Unit  Reynolds  number 

r  Radial  coordinate 

S  Surface  area 

s*  Directed  surface  area 

T,,T2  See  Eq.  (23) 

t  Time 

U  See  Eq.  (1) 

UDif  Turbulence  model  function  (see  Eq.  (12)) 

u  Longitudinal  velocity 

Vol  Volume 

v  Radial  velocity 

x  Longitudinal  coordinate 

y  Normal  coordinate 

Ymax  Turbulence  model  function  (see  Eq.  (11)) 

7  Ratio  of  specific  heats 

Ar  Radial  increment 

At  Time  increment 

Ax  Longitudinal  increment 

5  Vorticity-gradient  relaxation  parameter  (see  Eq.  (13)) 

e  Eddy  viscosity 
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€R  Relaxed  eddy  viscosity 

X  Second  coefficient  of  viscosity 

H  First  coefficient  of  viscosity 

£  Tangential  coordinate 

q  Density 

a  Normal  stress  including  pressure 

a  Normal  stress  excluding  pressure 

t  Tangential  stress 

0  See  Eq.  (28) 

u?  Vorticity 

SUBSCRIPTS 
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c 
D 
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in 
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out 
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oo 


Denotes  mesh  cell  centroid 
Denotes  reference  to  model  chord 
Denotes  reference  to  model  diameter 
Numerical  index 

Denotes  reference  to  inner  portion 
Numerical  index 

Denotes  reference  to  outer  portion 
Deontes  radial  direction 
Denotes  wall  conditions 
Denotes  longitudinal  direction 
Denotes  circumferential  direction 
Denotes  free-stream  conditions 
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